Quantum phases of electric dipole ensembles in atom chips 
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We present how a phase factor is generated when an electric dipole moves along a closed trajec- 
tory inside a magnetic field gradient. The similarity of this situation with charged particles in a 
magnetic field can be employed to simulate condensed matter models, such as the quantum Hall 
effect and chiral spin Hamiltonians, with ultra cold atoms integrated on atom chips. To illustrate 
this we consider a triangular configuration of a two dimensional optical lattice, where the chiral spin 
Hamiltonian Si ■ Bj x ffk can be generated between any three neighbours on a lattice yielding an 
experimentally implementable chiral ground state. 
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The interaction of the electromagnetic field with 
charged particles has unveiled a variety of geometrical 
and topological effects, such as the Aharonov-Bohm ef- 
fect Q, and plays a predominant role in the generation 
of collective phenomena, such as the quantum Hall effect 
and fractional statistics 0. The Aharonov-Bohm effect 
associates a quantum phase with the state of a charged 
particle. This phase is proportional to the enclosed mag- 
netic flux when the particle undergoes a looping trajec- 
tory and it is obtained even if the magnetic field is van- 
ishing in the immediate neighbourhood of the trajectory. 
This intriguing fact inspired a series of extensions, e.g. 
the reciprocal or the dual of the original Aharonov-Bohm 
effect H, ^ IE 111- Here, we shall consider another pos- 
sible generalisation that can be applied to the neutral 
atom technology of Bose-Einstein condensates and pro- 
vides a unique laboratory for the simulation of condensed 
matter models, e.g. chiral spin lattices or the quantum 
Hall effect, exhibiting advanced controllability and long 
coherence times. Such engineered systems can play a sig- 
nificant role in probing topological quantum simulations 
and error- free quantum information processing 0. 

With advances in cold atom technology it has been 
possible to simulate many solid state effects 0, 0, 0| 
by employing neutral atoms and techniques from quan- 
tum optics. Initially, it seemed that the neutrality of the 
atoms significantly restricts the range of phenomena this 
technology can simulate. Nevertheless, great effort has 
been made to develop neutral atom techniques, generat- 
ing suitable phase factors by laser radiation to simulate 
the behaviour of charged particles _Q, .H, J^] . These 
include proposals for the creation of the Hofstadter but- 
terfly 12'| and the lattice implementation of the quantum 
Hall effect 0|. These realisations arc limited in terms 
of their possible applications and it still remains to de- 
velop a technique for the efficient simulation of the effect 
of electromagnetic flelds on charged spin particles. 

Here, an alternative approach is presented, which al- 
lows us to go well beyond these applications. In the 
following we consider the electric dipole moment of the 



atoms in the presence of an appropriate electromagnetic 
held. The resulting effect is equivalent to a charged par- 
ticle moving in the presence of a magnetic fleld. This can 
simulate both the continuous case of electrons confined in 
a two dimensional plane, resulting in the quantum Hall 
effect as well as the case of discrete lattice sys- 

tems, which can result in the realisation of chiral states 
01 . The advantage of the present scheme is based on the 
simplicity of the required control procedure: utilising the 
atomic electric dipole as an additional degree of freedom 
provides adequate resources for probing even more com- 
plex structures, such as charged spin systems. 

Let us first see the general setup for implementing in- 
teractions between a particle and an external electro- 
magnetic field. Consider the case where a particle has 
a charge e, and an electric or a magnetic moment, given 
by de or dm respectively. The minimal coupling of the 
particle with the electromagnetic field is given by substi- 
tuting its momentum for 



p^p + eA + dyny- E + {de- V)A, 



(1) 



where E is the electric field and A is the magnetic vector 
potential. The second term in requires that the parti- 
cle is charged and gives rise to the well known Aharonov- 
Bohm effect (ij. While its topological character with 
the employment of a magnetic solenoid is of great con- 
ceptual value, and has been verified experimentally |l8l |. 
it has also given rise to a variety of applications in the 
solid state arena |^. The third term in |^ is the ori- 
gin of the Aharonov-Casher effect Q, which is recipro- 
cal to the Aharonov-Bohm effect. It involves the cir- 
culation of a magnetic dipole around a charged straight 
line and it has been experimentally verified Re- 
cent experiments have been performed that generalise the 
Aharonov-Casher effect thereby partly overcoming orig- 
inal technological complications I^J- The fourth term 
involves the coupling of the electric dipole moment to a 
differential of the vector potential and its consequences 
in cold atom technology will be the focus of this Letter. 
In fact, electric dipoles can give rise to a variety of dif- 



ferent phenomena that generate quantum phases when 
they undergo a cyche trajectory C. Starting from 
and by employing Stokes's theorem the phase factor con- 
tribution to the final state of the dipole is eventually 
given by 



I [V X (B X 4) + V X V(4 • A)] 
J s 



ds, (2) 



where 5 is a surface bordered by the cyclic path C and 
ds is its elementary area. The second term on the right 
hand side is the curl of the gradient oi d^ ■ A and can be 
taken to be zero, assuming continuity. The case where 
de- A is a multi-variable function, as might be produced 
when the dipole passes through a magnetic sheet, results 
in topological effects that have been studied in Hence, 
from (0), we obtain 



[{de ■ W)B - de{V ■ B)] - ds. 



(3) 



The second term on the right hand side is zero as can be 
seen from the Maxwell equations. Alternatively, a con- 
figuration with an infinite chain of magnetic monopoles, 
that alter the usual Maxwell equations, can lead to the 
generation of a topological phase dual to the Aharonov- 
Casher one 0, . For V • B = we finally obtain 



Is 



V)B • ds. 



(4) 



This equation can be viewed as the differential Aharonov- 
Bohm effect. By inspection of relation we see that a 
nontrivial phase can be produced if we generate an in- 
homogeneous magnetic field in the neighbourhood of the 
dipole. In particular, a non-zero gradient of the magnetic 
field component perpendicular to the surface S, varying 
in the direction of the dipole, ensures a non-zero phase 
factor. The dual effect with a magnetic dipole and a 
gradient of an electric field can also give a similar phase 
factor that originates from the third term in while 
a homogeneous electric field will just orient the electric 
dipole appropriately. As a realisation of I^J , we can take 
S to lie on the x-y plane and de to be perpendicular 
to the surface S. A non-zero phase, is produced if 
there is a non-vanishing gradient of the magnetic field 
along the z direction (see Figure ^a)). Alternatively, if 
de is along the surface plane, then a non-zero phase is 
produced if the z component of the magnetic field has a 
non- vanishing gradient along the direction of de as shown 
in Fig. n^b). Furthermore let us consider the case where 
the gradient of the magnetic field is spatially restricted, 
e.g. in a finite interval along the x direction. Then a tra- 
jectory of an electric dipole outside this area, S*, gives rise 
to a topological phase independent of the local details of 
the trajectory. 
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FIG. 1: The path circulation of the electric dipole in the in- 
homogeneous magnetic field. Fig. (a) depicts magnetic field 
lines suflicient to produce the appropriate non- vanishing gra- 
dient of along the z axis. Fig. (b) depicts only and 
how it varies along the direction of the dipole, de. 



These configurations of magnetic field gradients have to 
satisfy the Maxwell equations, V • i? = and V x. B — Q. 
Indeed, for dB^jdx ^ one should allow for an addi- 
tional variation of the x component of the magnetic field 
along the z direction with equal gradient. It can be veri- 
fied easily that this satisfies Maxwell equations and pro- 
vides the nontrivial phase factor given in Q). The same 
holds in the presence of a uniform electric field used to 
orient the electric dipole appropriately. 

Let us now consider the implementation of the phase 
factor with neutral atoms confined on a two dimen- 
sional plane. The electric dipole moment of an atom, 
which can be induced with the application of an external 
homogeneous electric field, is of the order \de\ ~ n'^eao, 
where n is the principal quantum number of the electron 
orbit, e is its charge and oq is the Bohr radius. It is even 
plausible to consider the ernployment of Rydberg atoms 
for which n is quite large |22l |. By applying an exter- 
nal macroscopic magnetic field we achieve gradients of 
the order dB^/dx « 10^-10^ Gauss/m. Much larger gra- 
dients of the magnetic field are feasible with atom chip 
technology due to the miniaturisation of their structure. 
For example, gradients of the order of 10® Gauss/m have 
already been reported at Caltech for the manipulation 
of ultra cold neutral atoms . This technology should 
prove very fruitful for the future of Bose-Einstein exper- 
iments where whole atomic ensembles can be mounted 
on atom chips ,24 Esl |2^ . With these at hand we can 
calculate that for atomic number n = I and for the cir- 
culation of an area of the order of the optical wavelength 
(~ 1/im) squared one can achieve phases of the order of 
« 8 > 27r. Hence, by varying the gradient of the mag- 
netic field or even the orientation of the electric dipole 
via an additional electric field it is possible to obtain an 
arbitrary value for the phase factor e"^. This provides a 
great degree of controllability when simulating charged 
particles in the presence of magnetic fields of arbitrar- 
ily large intensity. Next we shall see applications of this 
property. 

An important tool in the manipulation of atomic en- 
sembles is the employment of optical lattices. They can 
generate one, two or three dimensional structures of po- 
tential minima for the atoms. In particular, a one dimen- 
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sional optical lattice is sufRcient to confine the atomic 
ensemble in parallel planes above an atom chip, provid- 
ing two dimensional confinement, a necessary condition 
for the implementation of the proposed phase factor Q 
with atomic ensembles. For a sufficiently high flux of 
the effective magnetic field through an area S, given by 
e*B*S/h, and low two-dimensional densities of atoms, 
n2D, one can obtain 1271] . for example, the fractional fill- 
ing factor 1/ — n2D / {e* B* / h) = 1/2. This corresponds 
to the fractional quantum Hall effect [ifil ll^ which can 
be described by the m = 2 Laughlin wavefunction |2^ . 

Moreover, additional optical lattices can be applied 
to create a regular structure on this plane in order to 
simulate two dimensional topological spin effects. This 
is achieved by generating complex tunneling interac- 
tions along the planar lattice sites. In particular, con- 
sider two species of atoms, namely a = {T,J.}, that 
are superimposed with particular configurations of op- 
tical lattices. The evolution of the system, for atoms 
restricted in the lowest Bloch-band, is described by the 
Bose-Hubbard Hamiltonian that is comprised of tunnel- 
ing transitions of atoms between neighbouring sites of 
the lattice, V = -J^iai'^i^lcr^^i+icr + H.c), and col- 
lisional interactions between atoms in the same site, 

= ■^J^icrcr' U<y<y"^L'4a'°'i<y'<^i'^ ^^^^ COUpliugS, C/|| , 

Uii and [/||. By arranging for the tunneling couplings 
to be small with respect to the collisional ones, the sys- 
tem can be brought into the Mott insulator phase with 
only one atom per lattice site ll^. Hence, one can as- 
sume that the space of states of the system is spanned 
by the basis states | t)i = \n] = = 0) and 

I l)i = |nj = 0,n'l — 1), where nf is the number of 
atoms of species a in site i. 

In the rotated frame with respect to H'^^'> one can ex- 
pand the total Hamiltonian in terms of the small tun- 
neling interactions, obtaining an effective Hamiltonian 
that describes the tunneling of atoms without populating 
energetically unfavourable states that have two or more 
atoms per site. Subsequently, this Hamiltonian can be 
expressed in terms of the Pauli operators. Such expan- 
sion has provided a variety of two and three spin Hamil- 
tonians 0, |^ 0. By the additional employment of a 
magnetic field gradient and by considering the electric 
dipole of the atoms, it is possible to generate complex 
tunnehng couplings of the form J = e'^'^\J\ with 



We also focus initially on the case where the optical lat- 
tices generate a two dimensional structure of equilateral 
triangles. Such a non-bipartite structure is necessary in 
order to manifest the breaking of the symmetry under 
time reversal, T, in our model, eventually producing an 
effective Hamiltonian that is not invariant under com- 
plex conjugation of the tunneling couplings. Moreover, 
as the second order perturbation theory is manifestly T 
symmetric, we need to consider the third order. In that 
case the effective Hamiltonian becomes 



(5) 



with a = (cr^ , (T^ , (7^ ) and (...) denoting nearest neigh- 
bour sites. The couplings appearing in (O are given by 



^(4) ^ , 
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TT 



(4 • V)v4 • dx. 



Extraneous Zeeman terms in the z direction can be dealt 
with by Raman transitions that effectively create a com- 
pensating field of the form Bz = 2{J^ /U]] — /Uii)- 
The coupling constants, r^*) , can take various values. For 
example, one can choose U^i — > oo (e.g. by employing 
fcrmionic atoms) so that t*-^-* vanishes. For certain values 
of the collisional couplings one can also vary the tunnel- 
ing couplings such that t^^^ — 0. Moreover, by appro- 
priately tuning the phase of the tunneling couplings one 
can choose either t'^'^^ or r'-'*-' to vanish. 

Remarkably, with this physical proposal, a chiral three- 
spin interaction appears in which can be isolated, 
especially from the Zeeman terms that are predominant 
in equivalent solid state systems. This interaction term 
is also known in the literature as the chirality operator 
. It breaks the time reversal symmetry of the system 
as a consequence of the externally applied field. A chiral 
spin state is then a state for which the expectation value 
of the chirality operator has a nonzero value independent 
of the position of the plaquette (j, j, k). 



Here Xi and Xi+i denote the positions of the lattice sites 
connected by the tunneling coupling J. This is equivalent 
to the tunneling of a charged particle (e.g. electron) along 
a lattice in the presence of a magnetic field. 

In order to isolate the new effects generated by the 
consideration of complex tunneling couplings, we restrict 
ourselves to purely imaginary ones, i.e. = ±i\J^\. 



As a particular example, for studying the behaviour 
of the three spin interaction term, we take a hexagonal 
configuration of 19 spin (see Fig. I3a)). On this lattice 
we simulate Hamiltonian (|3J) for r'""^' = 0, t^^^ = r • r^"*) 
(t'*' > 0) and r'^'^^ — with the additional condition that 
the spins on the boundary experience a strong magnetic 
field oriented in the z-direction. A numerical simulation 
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creases due to frustration. As can be easily calculated 




for a single triangle, He has doubly degenerate ground 
states. This degeneracy is actually lifted by the chiral 
interaction giving a non-zero x for t > 0. Indeed, one 
finds that on a triangle the common ground state of the 
two interactions is |*) = ;^(| TTD+^^I TiT)+^^^l ITT)) 
which has non-vanishing chirality. This fact is in agree- 
ment with the persistence of chirality for large r as shown 
in Fig. Eld), which is the regime of the adopted pertur- 
bation theory. This experimentally feasible domain of 
couplings exhibits as a ground state a chiral spin state 
that may be possible to detect with the state of the art 
technology. 

This work was supported by the Royal Society. 



FIG. 2: (a) The hexagonal structure with 19 spins on the 
vertices and the 24 plaquettes. (b) The energy gap, A, from 
the ground state to the excited one as a function of the cou- 
phng r. (c) The chirality, x, as a function of the plaquette 
positioned on the plane for three different couplings r cor- 
responding to (d), where the chirality, Xi of each plaquette 
is given as a function of the coupling r and the particular 
plaquette M. 



has been performed to obtain the energy gap. A, above 
the ground state (see Fig. [2b)) as well as the chirality 
of the ground state, x = ■ o^j x <^k), for each triangular 
plaquette of neighbouring sites i, j and k (see Fig. |5Jc) 
and (d)). Indeed, Fig. [2Ib) shows that there is criticality 
behaviour for « —0.7 and « —0.1, where the 
energy gap becomes zero. In Fig. EJd) the chirality for 
different values of r and for different plaquettes on the 
plane numbered by M is displayed where a reconstruction 
of the chirality on the plane is given inj^c,!), (c,2) and 
(c,3) for the three distinctive areas. 

For T < the exchange interaction. He — 

J2{ij){'^f'^j + dominates the chiral interaction. 

He = J2{ijk) ' ^ ^k, forcing the spins to be aligned, 
pointing upwards in agreement to the boundary condi- 
tions. In this case x is almost zero (see Fig. |2Ic,3)). 
This holds also when we consider the Ising or a Zeeman 
interaction dominating He- For Tc-^ < r, the interaction 
He dominates giving a non-zero value to x- At Tcj^ and 
the chiral order parameter jumps indicating a quantum 
phase transition towards a chiral phase of the spin sys- 
tem. In particular, the plaquettes that have two sites on 
the boundary have small x &s both of the boundary spins 
tend to have parallel orientation. There is a small jump 
in X around as He induces a transition between en- 
ergetically favourable states from ferromagnetic (r < 0) 
to antiferromagnetic ones (r > 0), for which case x i^" 
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